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© Verfahren und Vorrichtung 2ur Verkehrslokalisierung in einem zellularen Mobilfunknetz 

® Der Anmeldungsgegenstand betrifft Verfahren und 
Vorrichtung zur Verkehrslokalisierung in einem zellularen 
Mobilfunknetz, bei dem mit Hilfe eines Storungsmodells, 
eines Kundenmodells und einer Erwartungsmaximie- 
rung, wodurch Beobachtungs- und Obergangswahr- 
scheinlichkeiten fur ein Hidden-Markov- Model I ermittelt 
werden, umdie Wahrscheinlichkert, mit der ein bestirnm- 
ter Kunde einen Bericht von einem bestimmten Pixel aus 
sendet, zu ermitteln. Ausgehend von diesem Hidden-Mar- 
kov-Modell konnen auch hochstwahrscheinliche Pfade ei- 
nes oder mehrerer Kunden ermittelt werden. 
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Beschreibung 

[0001] Die Erfindung betriffi ein Verfahren zur Verkehrslokalisierung in einem zellularen Mobilfunknetz, bei dcm 
mindestens eine Mobilstation, mindestens eine bedienende Basisstation und Nachbarbasisstauonen vorhanden sind, wo- 
5 bei mindestens eine Sequenz aufeinanderfolgender Berichte mit Empfangsstarken an eine Recheneinheit gemeldet wer- 
den und wobei ein zu bedienender Bereich in Rasterzellen unterteilt isL 
[0002] Ein derartiges Verfahren ist aus der intemationalen Anmeldung WO 98/15149 bekannt. 
[0003] Prinzipiell kann eine genaue Lokalisierung beispielsweise durch einen itn Mobilteil vorhandenen GPS-Emp- 
ranger relativ genau erfolgen. Dies ist jedoch teuer undbenougt eine freie Sicht zu drei Satellites Neben dem aufwen- 

10 digeren Mobilteil sind auch groBere Veranderungen in den Basisstationen erforderlich. 

[0004] Die intrazellulare Verkehrsverteilung in Mobilfunknetzen kann nur indirekt beobacfatet werden, indem man die 
von den Mobilslaiionen gemessenen RXLEV- und unter Umstanden TA-Werte (timing advance) sammelL Die RXLEV- 
Werte sind eine Anzahl von codierten Empfangsstarken der Basisstation und der nachsten Nachbarstationen, die fur ei- 
nen sogenannten Handover in Frage kommen. Die TA-Werte entsprechen der zeitJichen Verzogerung mit der die Bursts 

15 an der Basisstation empfangen werden und hangen von der Entfernung des Mobilteils von der Basisstation ab. Diese 
RXLEV- Werte und unter Umstanden auch die TA-Werte sind jedoch wegen Fading- und Shadowing-EfTekten stark ver- 
rauscht 

[0005] In den IEEE Transactions on Vehicular Technology, Vol. 48, No. 5, Sept. 1999, Seiten 1558-1562 sind einige 
Methoden, zum Beispiel Kalmanfilter, Neuronale Netzwerke oder Mustererkennungsverfahren auf der Basis von Hid- 
20 den-Markov-Modellen als prinzipielle Losungsansatze bei komplizierten Umgebungsbedingungen mit Shadowing usw. 
erwahnL 

[0006] Die nun zugrunde liegende Aufgabe besteht nun darin, ein Vsrfahren bzw. eine Vbrrichtung anzugeben, bei 
dem/der, auch bei starken Storungen durch Fading und Shadowing, eine schnelle und vor aDem eine genauere Lokalisie- 
rung der einzelnen Mobilstation und damit eine bessere und effizientere Analyse des Kundenverhaltens insgesamt mog- 
25 lichist 

[0007] Diese Aufgabe wird hinsichllicb des Verfahrens durch die Merkinale des Patenlanspruchs 1 oder des Patentan- 
spruchs 6 und hinsichtlich der Vorrichtung durch die Mcrkmalc des Anspruchs 10 oder 11 encicht 
[0008] Vorteilhafte Ausgestaltungen der Erfindung ergeben sich aus den weiteren Anspruchen. 
[0009] Ein besonderer Vorteil der Erfindung liegt in seiner hohen Aufldsung, die bei Simulationen zwischen 20 m und 
30 1 00 m lag. Die Erfindung eignet sich vor allem auch fur niedrige Geschwindigkeiten, wie sie beispielsweise bei FuBgan- 
gern auftreten, oder fur Bewegungen mit relativ geringer Distanz. 

[0010] Die Erfindung weist insbesondere beim Auffinden von sogenannten "hot spots" in der Verkehrsverteilung Vor- 
teile auf, da automatisch mehrere Kundenpfade uberlagert werden und eine Tendenz besteht, daB Kunden an typischen 
Kundenpfaden gesammelt werden. 

35 [0011] Die Grundidee der Erfindung besteht darin, eine Mobilfunkzelle in Rasterzellen zu unterteilen, ein Stormodell 
auf der Basis einer approximierte GauBverteilung mit Komelationen zu verwenden, um die Unterschiede, aufgrund von 
Fading und Shadowing, zwischen dem berichteten und vorhergesagten Feldstarken zu modellieren, ein Hidden-Markov- 
Modell fur die Rasterzellen anzugeben, das die Bewegungen der Mobilfunkteilnehmer bzw. der Kunden beschreibt und 
die Ubergangswahrscheinlichkeiten des Hidden-Markov-Modells dadurch zu ermitteln, dassder Erwartungswert der be- 

40 richteten Feldstarken maxirniert wird. Die stationare Verkehrsverteilung kann aus den tJbefgangswahrscheinlichkeiten 
des Hidden-Markov-Modells errechnet werden. 

[0012] Aus den Ubergangswahrscheinlichkeiten des Markov-Modells und den Beobachtungswahrscheinlichkeiten des 
Stormodells kann aus einer Sequenz von Berichten ein hochst wahrscheinlicher Pfad eines Mobiiteilnehmers ermittelt 
werden. 

45 [0013] Die Erfindung wird nun anhand eines in den Zeichnungen dargestellten Ausfuhmngsbeispiels naher erlautert. 
[0014] Dabeizeigt 

[0015] Fig. 1 eine Anordnung einer Mobilstation und mehreren Basisstationen zur Erlauterung der Erfindung und 
[0016] Fig. 2 einen Ausschnitt aus einer durch Rasterzellen unterteilten Mobilfunkzelle in Form von drei solchen Ra- 
sterzellen. 

50 [0017] In Fig. 1 ist eine Mobilstation MS, eine diese Mobilstation bedienende Basisstation BS sowie zwei zur Mobil- 
station nachst gelegene Basisstationen NS1 und NS2 dargestellt Die Mobilstation MS misst Feldstarkewerte RXLEV1 
der Basisstation BS Feldstarkewerte RXLEV2 der Station NS1 und Feldstarkewerte RXLEV3 der Station NS2 und mei- 
det diese gemessene Feldstarkewerte zugeordnet zu den jeweiligen Basisstationen in Form einer Meldung m" direkt an 
eine Recheneinheit RE oder altemativ in Form einer Meldung m' an eine Basisstation, z. B. die Basisstation BS, meldet, 

55 welche dann z. B. Meldungen sammelt, unter Umstanden weitere MeB werte hinzufugt und in Form einer Gesamtmel- 
dung m an die Recheneinheit RE meldet Diese Recheneinbeit RE wird sich typischerweise in einer Basisstation befin- 
den, kann aber auch in einer von den Basisstation raumlich getrennten Einheit oder sogar in einer mobilen Station vor- 
handen sein. Mindestens eine zu der bedienenden Basisstationen BS, gehorige Mobilfunkzelle ist durch ein Raster P mit 
einer Vielzahl von Raster oder Pixelzellen PIX weiter unterteilt. Die GroBe einer solchen Rasterzelle liegt dabei typi- 

60 scherweise zwischen 1 0 x 1 0 m und 50 x 50 m, wobei dies beispielsweise davon abhangt, ob es sich um einen stadtischen 
oder landlichen Bereich handelL Das Raster P bildet die Grundlage fur ein Hidden-Markov-Model, das in der Rechen- 
einheit RE implementiert isL 

[0018] In Fig. 2 sind drei solcher Rasterzellen PIX am Beispiel von drei benachbarten Rasterzellen pi, p2 und p3 dar- 
gestellt. Zwischen diesen Rasterzellen sind Obergangswahrscheinlichkeiten u vorhanden, wobei zwischen pi und p2 
65 eine Ubergangswahrscheinlichkeit ul2 und zwischen den Zellen p2 und p3 eine Obergangswahrscheinlichkeit u23 vor- 
handen ist. Falls die Zelle pi die erste Zelle eines vom Mobilfunkteilnehmer beschrittenen Pfades darstellt^ besteht cine 
Eingangswahrscheinlichkeit uOl zum Obergang auf die erste Zelle pi. In entsprechender Weise, wenn die Zelle p3 eine 
letzte Zelle eines vom Mobilfunkteilnehmer beschrittenen Pfades darsteilt, besteht fur das Vferlassen der Zelle p3 eine 
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Ausgangswahrscheinlichkeit ukO. Dariiber hinaus beslehen fur jedes Pixel Beobachtungswahrscheinlichkeiten f« die zum 
Ausdruck bringen, mil welcher Wahrscheinlichkeit ein bestimrnter Bericht von der jeweiligen Rasterzelle siammL 
[0019] Die Beobachtungswahrscheinlichkeit f p (m) = 0 fiir nicht in Fragekommende Pixel, weil z. B. der TA-Wert diese 
Pixel ausschlieBt, und berechnet sich sonst wie folgt 



fp(m) := i exp(- 



x 



-1 



x / 2 



wobei der Vektor x und der iransponierte Vektor x T aus den Differenzen von berichteten Feldstarkewerten aus einem je- 
weiligen Bericht m und zugehorigen vorausgesagten Feldstarkewerten gebildet werden. Die Matrix V" 1 ist das inverse 
einer Kovarianzmatrix V, die folgende Gestalt hat: 



a + b 



V = 




a + b 



wobei die Parameter a durch die physikalische Umgebung der jeweils betracbteten Basisstation und der Parameter b das 
phvsikalische Umfeld der Mobilslalion berucksichligL Die Konslante M dient lediglicb zur Nonnalisierung. Durch die 
ohcn ^cnanntcn Glcichung wird ein Storungsmodcll fur RastcrzcUcn gebildet, mit dcm fur Rastcrzcllcn allcr rclcvantcn 
Bcrichic Bcobachujngswahrscbeinlichkeiten f ermittelt werden. 

|(K)20| Das Verhalten eines typischen Kundens oder Benutzers einer Mobilstation wird als Markov- Kette modelliert 
Der /usiandsraum dieser Kette sind Rasterzellen oder Pixels eines Rasters P einer betracbteten Mobilfunkzelle und der 
spczicllc Zustand Off. Der Kunde kann sich von einem Off-Zustand zu einem Rasterpixel, der einen beginnenden Anruf 
darsiclli, bewegen. Der bei einem Pixel befindliche Kunde kann bei diesem Pixel verbleiben, sich zu einem Nachbarpixel 
bewegen oder zu einem Off-Zustand, der das Ende eines Telefongesprachs anzeigt, zuruckkehren. Normalerweise gibt es 
nur vicr NachbarpixeL Nord, Ost, Siid und West Falls sich die Mobilstation innerhalb einer Messzeit von bspw. 0,48 s 
weiier als die GroBe eines Rasterelements bewegt, mussen mehr Nachbarpixel berucksichtigt werden. Dieses Modell ba- 
sicrt auf der Idee, dass jede Rasterzelle zu einem typischen Kunden verhalten gehort. Zum Beispiel bewegen sich Mobil- 
siuiioncn innerhalb eines Buros oder Wohngebaudes ielativ selten, innerhalb FuBgangerbereichen bewegen sich Mobil- 
siaiioncn bzw. Kunden relativ langsam zu Nachbarpixeln und auf AutostraBen bewegen sich Mobilstationen haufig zu 
besiinimlen Nachbarpixeln in Abbangigkeit der Richtung der Strafie. 

1 002 1 J Ohne Vorwissen tiber das Kundenverhalten, abgesehen da von welche Pixel ausgehend von welchen Pixeln er- 
rcicht werden konnen, konnen aber zu Beginn Ubergangswahrscheinlichkeiten wie foigt angenommen werden: 

1 . u p>q : = l/(h + n), falls das Pixel q vom Pixel p wahrend einer Messperiode erreicht werden kann und sonst Up q = 
0. 

2. u Pw0 ff: = (h + n - n')/(h + n), fiir alle Pixel des Rasters, wobei n' die Zahl der Nachbarpixel von Pixel p ist. 
3- Uofrp = Up (0fi /Zp' e pu p ^o fiir alle Pixel des Rasters P und Uoff,off:=0 

n ist dabei die maximale Zahl von Pixeln, die von einem festgelegten Pixel aus erreicht werden kann und h kann bei- 
spielsweise 1 gewahlt werden. 

1 0022] Die l)bergangswahrscheinlichkeiten des Maitov-Modells konnen durch Training" aufgrund von Berichten der 
bcirachteten Mobilstationen verbessert bzw. den Verkehrsverhaltoissen im Mobilfunknetz angepasst werden. Dabei wird 
cin Erwartungswert fur die Berichte der beobachteten Kunden maximiert, das heiBt der groBte Wert fur folgenden Aus- 
druck gesucht: 



X 



number of cus tomer n path i length k 

n n v*. n fpj m 'j «w 

all possible p i% j customers t reports j 

[0023] Rine Variante des Baum-Welch-AlgorithmiK, der beispielsweise aus dem Buch voii Geoffrey J. McT^ughlan 
and Thriyambakam Krishnan mit dem Titel "The EM Algorithm and Extensions, Wiley series in probability and stati- 
stics, John Wiley and Sons, 1997, bekannt ist, lost dieses Erwartungswert-Maximierungswen-Problem auf besonders 
vorteilhafte Weise. Dieser Algorithmic versucht die Erwartung, dass der Kunde, welcher sich entsprechend der Markov- 
Keue verhalt, unter den spezielien Fading-Bedingungen diese Berichte mit den Messwerten erzeugt Dieser Algorithmus 
arbeitet iterativ. Er begihnt beispielsweise mit gleichfbrmig initialisierten Obei^angswahrscheinlichkeiten. Das Aus- 
gangsmodell korrespondiert mit der oben erwahnteo Situation, wo noch kein Wissen uber das Kundenverhalten vorhan- 
den ist Unter der Beschrankung, dass der Benutzer sich nur zu einem Nachbarpixel bewegen kann, berechnet der Algo- 
rithmus, welche tjbergange haufiger als andere sind, angezeigt durch Beobachtungswahrscheinlichkeiten f p (m) dass ein 
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Report m von einem besdmmten Pixel p slammt. Die nachste Iteration basiert auf den neuen verbesserten Obergangs- 
wanrscheinlichkeiten . 

[0024] Uber den BaunvWslch-Algorithmus werden einige Zwischenvariabien wie folgt definiert: 

>(*) ^ _ ..(*) 



15 




20 



/rWffl^ - (P) 



•«S'(p,») ~o 



30 



35 



a bzw, p sind dabei Wanxscheinlichkeiten, dass der i-te Benutzer Messwerte m^, . . bzw. m^, . . mi,! (ri) berich- 
tet und sich in Pixel p befindet, wenn er einen j-ten Bericht sendet, wobei l(ri) die Nummer des letzten Berichts in der i- 
ten Sequenz darsteilL Der Buchstabe 0 bedeutet den Off-Zustand. Die Variable % ist die Wahrscheinlichkeit, dass sich 
der i-te Benutzer sich beim Pixel p befindet wahrend er seinen j-ten Bericht sendet und sich damn zum Nachbarpixel q be- 
40 wegt Sind die Variablen £ berechnet, so sind die neuen Obergangswarirscheinbchkeiten wie folgt definiert 



45 



71 t=l 



50 



E, € j*w Eti J$S $<ft flO - 

* 

[0025] Es kann gezeigt werden, dass die neuen Obeigangswahrscheinlichkeiten eine erhohte Erwartung aufweisen. 
Durch Iteration der Berechnung der Zwischen variablen und das Updalen von Ubergangswahrscheinlichkeiten erhdht die 
60 Beobachtungswahrscheinlichkeiten. Es stellte sich heraus, das 10 bis 50 Tterationen fur eine angenaherte numerische 
Konvergenz iiblicherweise ausreichen. 

[0026] Die Zwischenvariabien £ in der letzten Iteration konnen aucb zur Berechnung der stationaren Verteilung der 
Markov- Kette ohne Mehraufwand herangezogen werden. Die Wahrscheinbchkeit Tip dass ein bestimmter Kunde einen 
Bericht von Pixel p sendet, ergibt sich aus folgender Formei: 

65 



4 
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—Steady— state 

[0027] Es kann nicht nur eine Wahrscheiniichkeit angegeben werden, mit der ein Kunde seinen Bericht von einem be- 
stimniten Pixel aus sendet, sondern es konnen auch hochst wahrscheinliche Pfade eines Kunden bzw. eines oder mehre- 
rer Mobilstationen ermittelt werden. Hierzu wird fur einen Pfad p b . . pt die Wahrscheiniichkeit, dass die Sequenz von 
Berichten mi, . . fQr ein festgelegtes Stdrungs- und Kundenmodell also fur feslgelegte Beobachtungswahrschein- 10 
iichkeiten und Ubergangswahrscheinlichkeiten maximiert wird Das Problem besteht also darin, einen Pfad p b . . ., p± zu 
finden, der die groBte Wahrscheiniichkeit hat, das heiBt den folgenden Wert maxirniert: 



• * 

path k*gthk 

H(Pi,p 2 ,...,p k ) = u^ ]\fAm j )u Pn 



P/ y J / Pj*Pj+\ 
reports j 



&j(p) ■= /pK)m{«p^i(?)} 
1>j(p) •= argm^c{«„ i9 <J i+l ( 9 )}. 

[0028] Der optimale Pfad p*i, . . ., p*i etgibt sicb durch ROckverfolgung der besten jeweiligen Auswahl durch 

p\ := argnm{tz fiiP 5i(p)} 

p* := tj-iip'.J fori = 2,...,/(r). 



15 



mit po: = pi+i: = 6 := off. Dieses Problem kann mit dem in der Nachrichtentechnik an sich wohl be kann ten Viterbi-Al- 20 
gorithmus gelost werden. Dieser Algorithmic wird hier wie folgt angewendet 

Eine Initialisierung fur eine gegebene Folge von Berichten r = (mi, , m*) erfolgt durch Setzen von 5k(p) = u p ^(m*) 

fur Pixel des Rasters. Dann erfolgt eine iterative Definition fur j = lr) - 1, . . ., und Pixel des Rasters. 



25 



30 



35 



[0029] Diese Informationen konnen genutzt werden, um die Qualitat der Dienste durch eine Optimierung der Einstel- 
lung zu verbessem, Problemzonen und -ursachen zu erkennen und zu beseitigen, notwendige Erweiterungen der Hard- 
ware zu planen und, als Fernziel, auf prekare Lastsituationen in Echtzeit zu reagieren. Dariiberhinaus konnen ortsabhan- 40 
gige Dienste fur einzelne Kunden angeboten werden. 



Patentanspriiche 

1. Verfahren zur Verkehrslokalisierung in einem zellularen Mobilfunknetz, 45 
bei dem mindestens eine Mobilstation (MS), mindestens eine bedienende Basisstation (BS) und Nachbarbasissta- 
tionen (NS1, NS2) vorhanden sind, 

mindestens eine Sequenz aufeinanderfolgender Berichte mit Empfangsstarken (RXLEV1 . . RXLEV3) an dem min- 
destens einen Mobilteil und/oder an der bedienenden Basisstation an eine Recheneinheit (RE) gemeldet werden, 
bei dem mindestens ein durch die mindestens eine Basisstation bedienter Bereich in einzelne Rasterzellen (PDC, pi, 50 
p2, p3) unterteilt wird, 

bei dem mit Hilfe eines Stdrungsmodells fur Rasterzellen aller relevanten Berichte Beobachtungswahrscheinlich- 
keiten (f) ermittelt werden und 

bei dem Eingangs-, Ausgangs- und/oder Ubergangswahrscheinlichkeiten fur die Rasterzellen der art bestimmt wer- 
den, dass ein Erwartungswert fur die mindestens eine Sequenz von aufemanderfolgenden Berichten durch veran- 55 
dern der Eingangs-, Ausgangs- und/oder Obergangs wahrscheinlichkeiten maximiert wird, wobei der zu maximie- 
render Erwartungswert dadurch definien ist, 

dass fur eine Sequenz aufeinanderfolgender Berichte und fur alle Pfade der selben Lange (k) das Produkt aus den 
Beobachtungswahrscheinlichkeiten (f) und den jeweiligen angenommenen Ubergangswahrscheinlichkeiten (u) 
zwischen den Rasterzellen des Pfades gebildet wird, 60 
das alle diese Produkte mit einer Eingangs- und Ausgangsubergangswahrscheinlichkeit (uOl , ukO) ftir die erste bzw. 
letzte Rasterzelle des Pfades zu einem erslen Gesamtprodukt multipliziert werden, 

dass die ersten Gesamtprodukte fur alk Sequenzen von aufeinanderfolgenden Berichten miteinander zu einem 
zweiten Gesamtprodukt multipliziert werden und 

dass die zweiten Gesamtprodukte fur alle moglichen Pfade zum jeweiligen Erwartungswert aufsummiert werden. 65 

2. Verfahren nach Anspruch 1 , bei dem der Erwartungswert fur die Berichte nach einer Variante des Baum-Welch- 
Algorithmus maximiert wird. 

3. Verfahren nach Anspruch 2, bei dem die Verkehrverteilung (rc) aus Zwischenwerten des Baum-Welch-Algorith- 



5 
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mus errechnet wird. 

4. Verfahren nach einem der Anspruche 1 oder 2, bei dem die Verkehrverteilung (rc) und/oder fur Rasterpunkte die 
ungefahre Geschwindigkeitsverteilung aus den Eingangs-, Ausgangs- und Ubergangswahrscheinlichkeiten (u) er- 
rechnet wird 

5. Verfahren nach einem der Anspriiche 1 oder 2, bei dem aus einer jeweiligen Distanz und der Zeit zwischen den 
Rasterelernenten zweier aufeinanderfolgender Berichte eine jeweilige ungefahre Geschwindigkeit der Mobilstation 
errechnet wird. 

6. Verfahren zur Verkehrslokalisierung in einem zellularen Mobilfunknetz, 

bei dem mindestens eine Mobilstation (MS), mindestens eine bedienende Basisstation (BS) und Nachbarbasissta- 
lionen (NS1, NS2) vorhanden sind, 

mindestens eine Sequenz aufeinanderfolgender Berichte nut Empfangsstarken (RXLEV1. . RXLEV3) an der min- 
destens einen Mobilstation und/oder an der bedienenden Basisstation an eine Recheneinheit (RE) gemeldet werden, 
bei dem mindestens ein durch die mindestens eine Basisstation bedi enter Bereich in einzelne Rasterzellen (PTX, pi, 
p2, p3) unterteilt wird, 

bei dem fur alle fur einen jeweiligen Bericht relevanten Rasterzellen BeobachmngswahischeinlichkeiteD (f) vorhan- 
den sind, 

bei dem fur alle fur den jeweiligen Bericht relevanten Rasterzellen Eingangs-, Ausgangs- und Ubergangswahr- 
scheinlichkeiten zwischen den Rasterzellen vorhanden sind und 

bei dem ein grdBtwahrscheiriUcher Pfad zu einer Sequenz aufeinanderfolgender Berichte dadurch gefunden wird, 
dass eine Gesamtwahrscheinlichkeit durch Verandern der Pfade maximiert wird, wobei die zu maximierende Ge- 
samtwahrscheinlichkeit eines Pfades dadurch definiert ist, 

dass fur alle aufeinanderfolgenden Berichte und fur alle Pfade der selben Lange (k) das Produkt aus den Beobacb- 
lungswahrscbeinlichkeiten (f) und den jeweiligen angenommenen Obergangswahrscheinlichkeiten (u) zwischen 
den Rasterzellen des Pfades gebildet wird und 

das alle diese Produkte zur Bildung der Gesamtwahrscheinlichkeit mit einer Eingangs- und Ausgangsubergangs- 
wuhrscheinlichkeil (uOl, ukO) fur die ersle bzw. letzte Rasterzelle des Pfades muldpliziert werden. 

7. Verfahren nach Anspruch 6, bei dem der grd£twahrschcinlichcn Pfad einer Sequenz von aufeinanderfolgenden 
Berichte der mindestens einen Mobilstation durch eine Van ante des Viterbi-Algorithmus gebildet wird 

8. Verfahren nach einem der Anspruche 6 oder 7, bei dem fur einzelne oder mindestens eine Gruppe von Sequenzen 
von aufeinanderfolgenden Berichten und Mobilstationen uber eine bestimmte Zeitdauer groBtwahrscheinliche 
Pfade der jeweiligen Sequenz und Mobilstation getrennt erf asst werden. 

( J. Verfahren nach einem der vorhergehenden Anspruche, bei dem mit Hilfe eines Storungsmodells fur eine jewei- 
lige in Frage kommende Rasterzelle eine Beobachtungswahrscheinlichkeit (f) derart ermittelt wird, 
dass Differenzen von gemessenen Empfangsstarken und vorhergesagten Empfangstarken fur eine jeweilige Raster- 
zelle zu einem Vektor (x) zusammengefasst werden, 

dass das Inverse einer Kovarianzmatrix, die die typischen zweiten Momente des Vektors (x) an der jeweiligen Ra- 
sterzelle charakterisiert, von links mit dem transponierten und rechis mit dem Vektor selbst sowie mit dem Faktor 
- 1 12 muldpliziert wird, 

dass fur das so erhaltene Produkt der zugehorige Wert einer Exponentialfunktion berechnet und dieses Ergebnis an- 
schlieBend durch eine Konstante normalisiert wird. 

10. Vorrichtung zur Verkehrslokalisierung in einem zellularen Mobilfunknetz, 

bei dem mindestens eine Mobilstation (MS), mindestens eine bedienende Basisstation (BS) und Nachbarbasissta- 
lionen (NS1, NS2) vorhanden, 

mindestens eine Sequenz aufeinanderfolgender Berichte mit Empfangsstarken (RXLEV1. . RXLEV3) an der min- 
destens einen Mobilstation und/oder an der bedienenden Basisstation an eine Recheneinheit (RE) gemeldet werden, 
bei der mindestens ein durch die mindestens eine Basisstation bedienter Bereich vorhanden ist, der in einzelne Ra- 
sterzellen (PIX, pi, p2, p3) unterteilt ist, 

bei dem ein Storungsmodells vorhanden ist, das fur alle Rasterzellen relevanter Berichte Beobachtungswahrschein- 
lichkeiten (f) ermittelt, und 

bei der ein Mittel zur Bestimmung von Eingangs-, Ausgangs- und/oder Ubergangswahrscheinlichkeiten zwischen 
den Rasterzellen derart vorhanden ist, dass ein Erwartungswert fur die Berichte der mindestens einen Mobilstation 
(MS) durch Verandern der Eingangs-, Ausgangs- und/oder Ubergangswahrscheinlichkeiten maximiert wird, wobei 
der jeweils zu maximierender Erwartungswert dadurch definiert ist, 

dass fur eine Sequenz aufeinanderfolgenden Berichte und fur alle Pfade der selben Lange (k) das Produkt aus den 
Beobachtungswahrscheinlichkeiten (f) und den jeweiligen angenommenen Obergangswahrscheinlichkeiten (u) 
zwischen den Rasterzellen des Pfades gebildet wird, 

das alle diese Produkte mit einer Eingangs- und Ausgangstibergangswahrscheinlichkeit (uOl, ukO) fur die erste bzw. 
letzte Rasterzelle des Pfades zu einem ersten Gesamtprodukt multipliziert werden, 

dass die ersten Gesamtprodukte fur alle Sequenzen von aufeinanderfolgenden Berichten miteinander zu einem 
zweiten Gesamtprodukt. multipliziert werden und 

dass die zweiten Gesamtprodukte fur alle moglichen Pfade zum jeweiligen Erwartungswert aufsummiert werden. 

11. Vorrichtung zur Verkehrslokalisierung in einem zellularen Mobilfunknetz, 

bei der mindestens eine Mobilstation (MS), eine bedienende Basisstation (BS) und mindestens eine Nachbarbasis- 
station (NS1,NS2) vorhanden sind, 

mindestens eine Sequenz aufeinanderfolgender Berichte mit Empfangsstarken (RXLEV1. . RXLEV3) an der rnin- 
destens einen Mobilstation und/oder an der bedienenden Basisstation an eine Recheneinheit (RE) gemeldet werden, 
bei der mindestens ein durch die mindestens eine Basisstation bedienter Bereich vorhanden ist, der in einzelne Ra- 
sterzellen (PIX, pi, p2, p3) unterteilt ist, 
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bei der fur alle fur einen jeweiligen Bericht relevanten Rasterzellen Beobachtungswahrscheiniichkeiten (f) vorhan- 
den sind, 

bei dem fur alle fur den jeweiligen Bericht relevanten Rasterzellen Engangs-, Ausgangs- und Ubergangswahr- 
scheinlichkeiten zwischen den Rasterzellen vorhanden sind und 

bei dern ein Mittel zum Finden eines groBtwahrseheinlichen Pfades der mindestens einen Mobilstation vorhanden 5 
ist, das den jeweiligen groBtwahrseheinlichen Pfad dadurch findet, dass eine Gesarntwahrscheinlichkeit durch \fer- 
andern der Pfade maximiert wird, wobei die jeweils zu maximierende Gesarntwahrscheinlichkeit dadurch definiert 
ist, 

dass fur aUe aufeinanderfolgenden Berichte und fur alle Pfade der selben Lange (k) das Produkt aus den Beobach- 
tungswabrscheinlichkeiten (f) und den jeweiligen angenommenen Ubergangswahrscbeinlichkeiten (u) zwischen to 
den Rasterzellen des Pfades gebildet wird und 

das alle diese Produkt zur Bildung der Gesamtwahrscheinlichkeit mit einer Eingangs- und Ausgangsiibergangs- 
wahrscheinlichkeit (uOl , ukO) fur die erste bzw. letzle Rasteiv^lle des Pfades multipliziert werden. 
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The invention relates to a method for traffic localization in 
a cellular mobile radio network, with at least one mobile 
5 station, one serving base station and adjacent stations being 
present, with at least one sequence of successive reports with 
reception strengths being signaled to a data processor and 
with an area to be served being subdivided into grid cells. 

10 A method of this kind is known from International Application 
WO 98/15149. 

In principle, a precise localization can, for example, take 
place relatively accurately by means of a GPS receiver in the 
15 mobile part. This is, however, expensive and requires a clear 
view to three satellites. In addition to the more expensive 
mobile part, considerable changes to the base stations are 
also necessary. 

20 The intracellular traffic distribution in mobile radio 

networks can be only indirectly observed, in that the RXLEV, 
and under certain circumstances TA (timing advance) , values 
measured by the mobile stations are collected. The RXLEV 
values are a number of coded reception strengths- of the base 

25 station and of the next adjacent stations that are used for a 
handover. The TA values correspond to the time delay with 
which the bursts are received at the base station and depend 
on the distance between the mobile part and the base station. 
These RXLEV values, and under certain circumstances also the 

3 0 TA values are, however, very noisy due to fading and shadowing 
effects. 

In IEEE Transactions on Vehicular Technology, Vol. 48, No. 5, 
Sept. 1999, pages 1558-1562, some methods, for example Kalman 
35 filters, neural networks or pattern detection methods based -on 
Hidden-Markov models are mentioned as the principle approaches 
to solving complicated environmental conditions with 
shadowing, etc . 



t 
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From the IEEE paper "New Methods for Processing GSM Radio 
Measurement Data: Applications for Locating, Handover and 
Network Management", Vehicular Technology Conference, 1994 
5 IEEE 44th Stockholm, SE, 8-10.06.1994; M. Junius et al . , 
methods are known from traffic localization in a cellular 
mobile radio network, whereby the cells are divided in sectors 
with a route section being described, for example, by a 
Hidden-Markov model, with characteristic changes in the 
10 reception strengths being represented by the parameters of 

this model and with reception strengths generated for certain 
route sections being compared with measured reception 
strengths and from this a route section more likely used by 
mobile subscribers being determined. 

15 

From the German application paper laid open for public 
inspection DE 199 34 845 Al , a method is known for estimating 
Hidden-Markov models with a low memory requirement. 

20 The object of the invention is now to provide a method or 

device whereby a fast and above all more accurate localization 
of the individual mobile station and thus a better and more 
efficient analysis of the customer behavior overall is 
possible even where there is strong interference due to fading 

2 5 and shadowing. 

The objective is achieved with regard to the method by the 
features of Claim 1 and with regard to the device by the 
features of Claim 10. 

30 

Advantageous embodiments of the invention are provided by the 
further claims. 

A particular advantage of the invention is in its high 

3 5 resolution, that in simulations was shown to be between 2 0 m 

and 100 m. The invention is also particularly suitable for low 
speeds, such as occur with pedestrians, or for movement over a 
relatively short distance. 



The invention has advantages, particularly when finding hot 
spots in the traffic distribution, because several customer 
paths are automatically superimposed and there is a tendency 
for customers to be collected on typical customer paths. 

The basic idea of the invention is to divide a mobile radio 
cell into grid cells, to use an interference model based on an 
approximated Gaussian distribution with correlations, in order 
to model the differences between the reported and predicted 
field strengths on the basis of fading and shadowing, to 
provide a Hidden-Markov model for the grid cells that 
describes the movements of the mobile radio user or of the 
customers and in this way to determine the transition 
probabilities of the Hidden-Markov model so that the expected 
value of the reported field strengths is maximized. The 
stationary traffic distribution can be calculated from the 
transition probabilities of the Hidden-Markov model. 

From the transition probabilities of the Markov model and the 
observed probabilities of the interference model, a maximum 
probability path of a mobile user can be determined from a 
sequence of reports. 

The invention is now explained in more detail by means of an 
exemplary embodiment shown in the drawings . 

These are as follows: 

Figure 1 An arrangement of a mobile station and several base 

stations for explaining the invention, and 

Figure 2 A section from a mobile radio cell subdivided by grid 

cells, showing three such grid cells. 

Figure 1 shows a mobile station MS, a base station BS serving 
this mobile station and two base stations NS1 and NS2 next to 
the mobile station. The mobile station MS measures the field 
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strength values RXLEV1 of the base station BS, field strength 
values RXLEV2 of the station NS1 and field strength values 
RXLEV3 of station NS2 and signals these measured field 
strength values, allocated to the relevant base stations, in 
5 the form of a message m" " directly to a data processor RE, or 
alternatively in the form of a message irT to a base station, 
e.g. base station BS, that then, for example, assembles 
messages, and under certain circumstances adds further 
measured values, and signals these to the data processor RE as 

10 a total message m. This data processor RE is typically located 
in a base station but can also be a unit spatially separated 
from the base station or can even be installed in a mobile 
station. At least one of the mobile radio cells belonging to 
the base station BS to be served is further subdivided by a 

15 grid P with a number of grids or pixel cells PIX. The size of 
such a grid cell in this case is typically between 10x10 m and 
50x50 m, whereby this, for example, depends on whether it is 
an urban or country area. Grid P forms the basis for a Hidden- 
Markov model implemented in the data processor RE. 

20 

Figure 2 shows three such grid cells PIX using three adjacent 
grid cells pi, p2 and p3 as an example. Transition 
probabilities u are present between these grid cells with a 
transition probability ul2 between pi and p2 and a transition 

25 probability u23 between cells p2 and p3 . If the cell pi is the 
first cell of a path used by a mobile radio user, there is an 
entry probability uOl of transition to the first cell pi. In a 
corresponding manner, if cell p3 is the last cell of a path 
used by a mobile radio user, there is an exit probability ukO 

30 of leaving cell p3 . Furthermore, there are observed 

probabilities f for each pixel, that express the probability 
of a specific report originating from the relevant grid cell. 

The observed probability f p (m) = 0 for pixels not under 
35 consideration, because, for example, the TA value precludes 
these pixels, is otherwise calculated as follows 

fp(m) := — exp ( -x T • V" 1 • x/2) , 

M 
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with the vector x and the transposed vector x T being formed 
from the difference between reported field strength values 
from a particular report m and the associated predicted field 
5 strength values. The matrix V" 1 is the inverse of a covariant 
matrix V, that has the following form 



with the parameters a being allowed for by the physical 
10 environment of the particular base station under consideration 
and parameter b by the physical environment of the mobile 
station. The constant M serves only for normalization. By 
means of the above equation, an interference model for grid 
cells is formed, with observed probabilities f being 
15 determined for all reports relevant to grid cells. 

The behavior of a typical customer or user of a mobile station 
is modeled as a Markov chain. The status space of this chain 
is grid cells or pixels of a grid P of a mobile radio cell 

2 0 under consideration and the special status OFF. The customer 

can move from an OFF status to a grid pixel that represents 
the start of a call. The customer at a pixel can remain at 
this pixel, move to an adjacent pixel or return to an OFF 
state that represents the end of a telephone call. Normally, 
25 there are only four adjacent pixels: North, East, South, and 
West. If the mobile station moves further than the size of a 
grid element within a measured period of, for example, 0.48 s, 
more adjacent pixels must be taken into account. This model is 
based on the idea that each grid cell belongs to a typical 

3 0 customer behavior. For example, mobile stations move 

relatively seldom within an office or dwelling, within 




6 

pedestrian areas mobile stations or customers move relatively 
slowly to adjacent pixels and on highways mobile stations 
frequently move to certain adjacent pixels depending on the 
direction of the road. 

5 

Without advance knowledge of the customer's behavior, 
regardless of which pixel can be reached starting from which 
pixel, transition probabilities can, however, be assumed from 
the start as follows: 

10 

1. u p , q : =1/ (h+n) , if the pixel q can be reached from pixel p 
during a measuring period and otherwise u p , q =0. 

2. u Pi0 ff : = (h+n-n" ) / (h+n) , for all pixels in the grid, with n" 
15 being the number of adjacent pixels to pixel p. 

3. u 0 ff, P = Up, off /Z P# e'pUp<,o for all pixels in grid p and 
u 0 f f ,of f : =0 . 

2 0 n in this case is the maximum number of pixels that can be 
reached from a specified pixel and h can, for example, be 
chosen as 1 . 

The transition probabilities of the Markov model can be 
25 improved, or matched to the traffic conditions in the mobile 
radio network, by "training" based on reports from the 
observed mobile stations. In this way, an expected value for 
the reports of the observed customers is maximized, i.e. the 
largest value for the following expression is sought. 

30 

number of cus tomer n path i leng th k 

x n *w, n fp,s m ^ u P,,^ 

all possible p t j customers i reports j 

A variant of the Baum-Welch algorithm known, for example, from 
the book by Geoffrey J. McLaughlan and Thriyambakam Krishnan, 
35 entitled "The EM Algorithm and Extensions, Wiley series in 
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probability and statistics, John Wiley and Sons, 1997, solves 
this expected value maximization value problem in a 
particularly advantageous manner. This algorithm attempts to 
anticipate that the customer who behaves according to the 
5 Markov chain generates these reports with the measured values 
under these special conditions . This algorithm works 
iteratively. It begins, for example, with the uniformly 
initialized transition probabilities. The initial model 
corresponds to the aforementioned situation where there is 

10 still no knowledge of the customer behavior. Under the 

limitation that the user can move only to the next adjacent 
pixel, the algorithm calculates which transitions are more 
frequent than others, indicated by observed probabilities f p (m) 
that a report m originates from a certain pixel -p. The next 

15 iteration is based on the new improved transition 
probabilities . 

By means of the Baum-Welch algorithm, certain intermediate 
variables are defined as follows: 



20 



c$(p) = «J2/,(r<,i), 




:= 0 

«W(r,)lP«»J — ~ (k) , .s (k) 

«Slr,)fei«) 0- 
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a or P in this case are probabilities that the ith user reports 
measured values mi,i f . . . , tru, j or mi, j, . . . , mi f i< r i) and is in pixel p 
when he sends a jth report, with l(ri) being the number of the 
last report in the ith sequence. The letter 0 means the OFF 
state. The variable £ is the probability that the ith user is 
at pixel p when he is transmitting his jth report and that he 
then moves to the adjacent pixel q. If the variables £ are 
calculated, the new transition probabilities are defined as 
follows . 



effort forpeP 



tx^ +1 > = n ifr) (?) forp€-PandgePU{0}. 

It can be shown that the new transition probabilities have an 
increased expectation. By iteration of the calculation of the 
intermediate variables and updating transition probabilities, 
15 the observed probabilities increase. It was shown that 10 to 
50 iterations were normally sufficient for an approximate 
numerical convergence . 

The intermediate variables £ in the last iteration can also be 
20 used without additional cost to calculate the stationary 

distribution of the Markov chain. The probability n p that a 
specific customer transmits a report from pixel p is obtained 
from the following formula. 

steady— state 

E£=i Ep'ep £ 9 epu{0} ft,i(p', q) 



25 



It is possible not only to give the probability with which a 
customer transmits his report from a certain pixel, but the 
most probable paths of a customer, or of one or more mobile 
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stations, can also be determined. Furthermore, for a path 
Pir • - . /Pk the probability that the sequence of reports mi, ... ,mk 
for a specified interference and customer model, i.e. for 
specified observed probabilities and transition probabilities, 
5 can be maximized. The problem therefore is to find a path 
Pi/.../Pk that has the greatest probability, i.e. that 
maximizes the following value: 

pat h leng th k 

H( Pl ,p 2 — /Pk )= u PotPi n/p, {m 3) u Pj .p j+1 

reports j 

10 

where p 0 : =Pk+i : =© : = off. 

This problem can be solved by the Viterbi algorithm, well 
knovm: in telecommunication engineering. This algorithm is in 
this case used as follows. 

15 

An initialization for a given sequence of reports r= (mi , . . . , m k ) 
takes place by setting 5 k (p) = u P/Q ff f P (m k ) for pixels in the 
grid. There then follows an iterative definition for j=l(r)- 
1, . . . , and pixels in the grid. 

20 

8j (p) := /p(mj) max{u pq S J+l (q)\ 
y/j (p) : = arg max }u p q 6 j + 1 (q)} . 

The optimum path p*i, ...,p* k is obtained by tracing the best 
25 selection in each case by means of 

Pi : = arg max {u 9/P 5 1 (p)} 

peP 

p* := ^j-itej-i) for j = 2, . . . ,/(r) . 

3 0 This information can be used to improve the quality of 

services by optimizing the setting, to detect and rectify 
problem zones and causes, to plan the necessary upgrades of 
hardware and, as a long-term aim, to react in real time to 
precarious load situations. Furthermore, local services can be 

35 offered for individual customers. 
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Method for traffic localization in a cellular mobile radio 
network, 

r- whereby at least one mobile station (MS) , at least one 
serving base station (BS) and adjacent stations (NS1, NS2) 
are present, 

- at least one sequence of successive reports with 
reception strengths (RXLEV1.. RXLEV3) is transmitted to a 
data processor (RE) at the mobile part, of which there is 
at least one, and/or at the serving base station, 

- whereby at least one area served by at least one base 
station is subdivided into single grid cells (PIX, pi, p2, 
p3), 

- whereby observed probabilities (f) for grid cells of all 
relevant reports are determined with the aid of an 
interference model and 

- whereby the entrance, exit and/or transition 
probabilities for the grid cells are determined in such a 
way that an expected value for at least one sequence of 
successive reports is maximized by changing the entry, exit 
and/or transition probabilities, with the expected value to 
be maximized being thus defined, 

for a sequence of successive reports and for all paths 
of the same length (k) the product of the observed 
probabilities (f) and the particular assumed transition 
probabilities (u) between the grid cells of the path is 
formed, 

all these products are multiplied by an entry and exit 
transition probability (uOl, ukO) for the first or last 
grid cell of the path to form a first total product, 

the first total products for all sequences of successive 
reports are multiplied by each other to form a second total 
product and 

the second total products for all possible paths are 
totalized to form the relevant expected value. 

Method in accordance with claim 1, 



whereby the expected value is maximized for the reports in 
accordance with a variant of the Baum-Welch algorithm. 

Method in accordance with claim 2, 

whereby the traffic distribution (71) is calculated from the 
intermediate values of the Baum-Welch algorithm. 

Method in accordance with one of claims 1 or 2 , 
whereby the traffic distribution (71) and/or for grid points 
the approximate speed distribution is calculated from the 
entry, exit and transition probabilities (u) . 

Method in accordance with one of claims 1 or 2 , 
whereby a particular approximate speed of the mobile 
station is calculated from a relevant distance and the time 
between the grid elements of two successive reports. 

Method for traffic localization in a cellular mobile 
network in accordance with claim 1, 

- with a path of greatest probability for a sequence of 
successive reports being found so that a total probability 
is maximized by changing the path, whereby the total 
probability of maximization of a path is defined in that, 
-- for all successive reports and for all paths of the same 
length (k) the product of the observed probabilities (f) 
and the relevant assumed transition probabilities (u) 
between the grid cells of a path is formed and 

all these products are multiplied by an entry and exit 
transition probability (uOl, ukO) for the first or last 
grid cell of the path in order to form the total 
probability. 

Method in accordance with claim 6, 

with the path of greatest probability of a sequence of 
successive reports of at least one mobile station being 
formed by using a variant of the Viterbi algorithm. 

Method in accordance with one of claims 6 or 7, 
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whereby for single or at least one group of sequences of 
successive reports and mobile stations the paths of 
greatest probability of the relevant sequence and mobile 
station are found separately over a specific time period. 

5 

9. Method in accordance with one of the preceding claims, 
whereby with the aid of an interference model for a 
relevant grid cell under consideration an observed 
probability (f) is determined in such a way 

10 that differences in measured reception strengths and 

predicted reception strengths for a particular grid cell 
are combined to form a vector (x) , 
-- that the inverse of a covariant matrix that 
characterizes the typical second moments of the vector (x) 

15 at the particular grid cell is multiplied from the left by 

* 

the transposed vector and from the right by the vector 
itself and by the factor -1/2, 

-- that the associated value of an exponential function is 
calculated for the product obtained in this way and the 
20 result is then normalized by a constant. 

10. Device for traffic localization in a cellular mobile radio 
network, 

- whereby at least one mobile station (MS) , at least one 
25 serving base station (BS) and adjacent base stations (NS1, 

NS2) are present in such a way that a data processing unit 
(RE) is present, to which at least one sequence of 
successive reports containing reception strengths (RXLEV1.. 
RXLEV3 ) is signaled at the mobile station, of which there 
3 0 is at least one, and/or at the serving base station, 

- whereby at least one area served by at least one base 
station is present, said area being subdivided into 
individual grid cells (PIX, pi, p2 , p3), 

- whereby an interference model for determining observed 

35 probabilities (f) for all grid cells of relevant reports is 

present , 
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- whereby a means for determining entry, exit and/or 
transition probabilities is present between the grid cells 
and 

- whereby means for maximizing an expected value for the 
reports of at least one mobile station (MS) are present, 
whereby this takes place by changing the entry, exit and/or 
transition probabilities and whereby the expected value to 
be maximized in each case is defined in that, 

-- for a sequence of successive reports and for all paths 
of the same length (k) the product of the observed 
probabilities (f) and the relevant assumed transition 
probabilities (u) between the grid cells of the path is 
formed, 

- all these products are multiplied by an entry and exit 
transition probability (uOl, ukO) for the first or last 
grid cell of the path in order to form a first total 
product , 

the first total products for all sequences of successive 
reports are multiplied by each other to form a second total 
product and 

the second total products for all possible paths are 
totalized to form the relevant expected value. 

1. Device for traffic localizing in a cellular mobile radio 
network in accordance with claim 10, 

- whereby a means for finding a path of greatest 
probability of at least one mobile station is present, that 
finds the particular path of greatest probability so that a 
total probability is maximized by changing the paths, 
whereby 

the particular total probability to be maximized is defined 
in that 

for all successive reports and for all paths of the same 
length (k) the product of the observed probabilities (f) 
and the particular assumed transition probabilities (u) 
between the grid cells of the path is formed and 
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--to form the total probability, all these products are 
multiplied by an entry and exit transition probability 
(uOl, ukO) for the first or last grid cell of the path. 



